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We formulate and apply a theory of electron-phonon interactions for the surface state of a strong 
topological insulator. Phonons are modelled using an isotropic elastic continuum theory with stress- 
free boundary conditions and interact with the Dirac surface fermions via the deformation potential. 
We discuss the temperature dependence of the quasi-particle lifetime in photoemission and of the 
surface resistivity. 
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Introduction. — One of the presently most active areas 
in physics is concerned with strong topological insulator 
(TI) materials [TJ [2] . In these systems strong spin-orbit 
couplings cause band inversion and a nontrivial topology 
of the map from momentum to Hilbert space 0] . In 
a TI, as long as time-reversal invariance remains unbro- 
ken, an odd number of massless surface Dirac fermion 
modes is guaranteed despite of the presence of a bulk 
gap A;,. The existence of metallic two-dimensional (2D) 
Dirac surface states has been convincingly established us- 
ing angle-resolved photoemission spectroscopy (ARPES) 
in bismuth selenides [TJ [5] . Typical reference materials 
are Bi 2 Te 3 or Bi 2 Se3, where ~ 0.3 eV allows to ob- 
serve these phenomena even at room temperature, but 
other material classes have also been predicted to possess 
a TI phase JT] . Attempts to observe electronic transport 
signatures of the surface state were only partially suc- 
cessful [BHH] since surface effects are easily masked by 
defect-induced residual bulk charge carriers. Anticipat- 
ing progress in achieving better purity, it is important to 
understand what intrinsically limits the surface conduc- 
tivity and the integrity of surface quasi-particles. Noting 
that the large and anisotropic static dielectric constant 
(e « 50 to 200 |9J implies a drastic reduction of direct 
Coulomb forces or charged impurity potentials, we here 
analyze consequences of the electron-phonon coupling on 
the TI surface state. 

We focus on long-wavelength acoustic phonons which 
dominate the physics at low energy scales. Previous work 
on the thermoelectric properties of Bi2Te3 has demon- 
strated that despite of the quintuple-layer crystal struc- 
ture, bulk acoustic phonons are reasonably well described 
as isotropic elastic continuum [TU1 HJJ, where the two 
Lame parameters of the theory determine the longitu- 
dinal and transverse sound velocities, c/ ~ 2800 m/s 
and Ct ~ 1600 m/s, respectively. Low-temperature elec- 
tronic transport is then limited by the deformation po- 
tential coupling to acoustic phonons, since piezoelectric 
couplings are suppressed by inversion symmetry [llj . TI 
experiments have so far only addressed the coupling to 
optical phonons [HI H2], cf. also studies for Bi surfaces 
[13] . However, massless 2D Dirac fermions are realized in 
graphene monolayers as well, where both theory |141 115] 
and experiment [16] have reported consistent results for 
the temperature (T) dependence of the resistivity (p): 



for T <C Tbg (with the Bloch-Griineisen temperature 
Tbg = 'ZhkpCs/kB, Fermi momentum hp, and sound ve- 
locity c s ), a p ~ T scaling is found, while p ~ T for 
T ^S> Tbg- Note that for a 2D electron gas with parabolic 
dispersion and dominant deformation potential coupling, 
one expects p ~ T 7 for T <C T B q [T7]. Below, we shall 
discuss the p(T) dependence of the TI surface state due 
to the coupling to acoustic phonons in detail. In addition, 
phonons are expected to cause quasi-particle decay. This 
implies, e.g., a finite ARPES linewidth [TS|[JJJ], where an 
anomalous behavior was observed in Bi2Se3 [2~0~j. 

In this paper, we formulate and study an analytically 
tractable effective low-energy theory for the TI surface 
state coupled to acoustic phonons. The phonon modes 
are obtained from isotropic elastic continuum theory in a 
half-space with stress-free boundary conditions j2TJ [22] . 
Their coupling to the surface fermions is predominantly 
via the deformation potential |23j . For concrete num- 
bers, we use published [TU1 [TTJ |2I] values for Bi 2 Te 3 . We 
compute the quasi-particle decay rate T and find T ~ T 
at high temperatures. This prediction should be observ- 
able by ARPES. Phonons also affect the surface resistiv- 
ity and yield a characteristic T-dependent contribution 
to p. Our theory is flexible enough to allow for the use 
of microscopically derived phonon modes and coupling 
matrix elements, e.g., resulting from (future) numerical 
force-constant calculations. 

Model. — We consider energies below the TI bulk gap 
where only surface electronic states are relevant. For the 
half-space z > 0, the surface-state wave function x( z ) 
follows from the low-energy band structure with Dirichlet 
boundary conditions at z — [25], 

-e— ), , ± = W|TMM 

(1) 

with normalization f\f and material parameters 
(So, -Mo, Mi) specified in Ref. [21]. Since M M X < 
and Bo /Mi > 0, we have Re(r/±) > and the state 
decays exponentially. One arrives at a massless 2D 
Dirac Hamiltonian (we set H = 1) [JJ, 

H e = € k s c{ s c ks , e ks = sv F \k\ - p, (2) 

fe,s=± 

with the Fermi velocity vf — 4.36 x 10 5 m/s and the 
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chemical potential p defining kp = \p\/vp. A helical 
eigenstate with helicity s = ± has its spin structure tied 
to the surface momentum k = (k xi k y ). Helical fermions, 
Cfe = (cfc + ,Cfc_) T , are connected to the usual spinful op- 
erators, dfc = (dfc^,d/4) T , by a unitary transformation, 



With K[ = \J\ — 7^ 2 and K t = \J\ — £ 2 , we obtain 
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where tan 9k = k y /k x . 

In order to describe noninteracting acoustic phonons 
we employ isotropic elastic continuum theory with stress- 
free boundary conditions at z = 0. We briefly summa- 
rize the resulting eigenmodes [211 E2] before turning to 
the electron-phonon coupling. Following the notation in 
Ref. |22j we label the modes by the quantum numbers 
A = (g,f2, A), with surface momentum q = (q x ,Qy)) fre- 
quency fl > 0, and mode type A € (H,T, L, R) explained 
below. In this non-standard but very convenient nota- 
tion, the frequency Q = Q\ is not specified in terms of q 
and A but represents a free parameter. With r = (x, y) 
and surface area A, the displacement field operator takes 
the form 

U(r, z,t)=J2 fo^o ^y^^ + h ' c " ^ 

y / 2p M Ail 

where b\ is a bosonic annihilation operator and pM — 
7860 kg/m 3 [10]. The noninteracting phonon Hamilto- 

man is H p = Ea^a( 6 W + V 2 )- Th e orthonormal 
eigenmodes ua(z) describe linear combinations oie ±lkl - tZ 



Using the above values for q.t we find £ ~ 0.92, ki ~ 0.85, 
Kt 0.39 and C ~ 1.20. 

Electron-phonon coupling. — The deformation poten- 
tial couples the local electron density to V • U(r,z), 
with a coupling constant a. Ref. [TT] gives the estimate 
a ~ 35 eV. This yields the second-quantized interaction 
Hamiltonian 
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where Uk in Eq. Q defines the matrix Xk q = U k+q U\. 
For the Rayleigh mode, the sum over Q should be omitted 
with the replacement O = c^q. With <f>^ specified in 
Eqs. |5| and ([6]), we obtain the e g -independent electron- 
phonon coupling matrix elements 
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waves, where k^t = \l l/ci^y — q 
zontal shear mode, A = H, with \\ e z x e q 
e q = q/q) decouples from all other modes and does not 
generate a deformation potential; hence it is not dis- 
cussed further. The remaining modes are given by 



with the electronic surface state \ m Eq. Q. For q <C 
Re (??_), the overlap integral above reduces to 4>^ [z = 0). 
First, the hori- i n wna t follows, we discuss physical consequences ob- 
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The incoming longitudinal mode, A = L, with a/ = 1 
and a t = 0, exists for f2 > c;q with real fc; t > 0. The 

has bi = — A and b t — B, where 



eigenstate <j>\ L t 



(where tained from the Hamiltonian H = H e + H p + H ep . In 
the concrete examples below, the chemical potential is 
p = vpkp = 0.05 eV, corresponding to the BG tempera- 
ture Tbg = 2kpCn/kB = 3.9 K and the Fermi tempera- 
ture Tp = 580 K. 

Lifetime broadening. — We begin with the self-energy 
E s (fe,w) for a helical eigenstate s = ±. Following stan- 
dard arguments [TH [TU], the main contribution is cap- 
tured to lowest nontrivial order in H ep . Noting that the 
"tadpole" diagram vanishes identically, the "rainbow" di- 
agram |18] yields independent contributions from each 
mode A, 
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The incoming transverse mode, A = T, with ai = and 
a t = 1, exists for Q > c t q. The eigenstate (f>\ fy has 
bi = —B and b t = —A. (For ctq < O < ciq, we have 
ki = i\ki\.) Finally, the energetically lowest solution is 
the Rayleigh surface wave, A = R, where a; = at = 
and ki t t — inuq. Here the dispersion relation is linear, 
O = cnq with surface velocity cr = fc f , i.e., f2 is not a 
free parameter in A anymore. Putting 7 = (c t /c;) 2 , we 
find 
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where n B (rip) is the Bose (Fermi) function. For the R 
mode, the g-integral has to include the additional factor 
<5(il — cnq), while for A = L,T, we have the respective 
constraint q < 0/c;, t . The decay rate T (A) = -2 Im E (A) 
describing lifetime broadening is then given by 
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FIG. 1. (Color online) Low-frequency behavior of the 



Eliashberg functions Fj,^ (SI) for the three relevant acous- 



tic phonon modes (k = 0.1/cf). In the rescaled units used 
here, the functions are approximately fc-independent. Inset: 
Same for the "transport" Eliashberg function T (see main 
text). 




FIG. 2. (Color online) Main panel: T-dependence of the 
decay rate V for k = kp and for k — 0.5fc_F- For k — 0.5fcF, 
only the L mode gives significant contributions. Inset: k- 
dependence of r for T = 3.9 K (solid line) and for T = 392 K 
(dashed line; the shown result has to be multiplied by 10). 



Here the Eliashberg function [18] is defined as 
<t ) (^)=E / ^\M$X kq>ss ,\ 2 6(u±n-e k+qtS ,) 

s' 

(11) 

which represents a phonon density of states weighted by 
the coupling matrix elements. Performing the angular 
integration yields the result 



" ks.oj 



(O) 



q± 



l 



2~nvpk 



qdq 



M. 



(A) 
qQ 



1+ 



q 2 - q 2 _ 



s'/2 



\fl + LU + vSl\ 



± k 



v F 



(12) 



where s' = s sgn(/i + u) + vSi). 

For a discussion of the lifetime, we now consider the 
on-shell case, ui — e ks . For the Rayleigh mode with cr -C 
Vp, we find s' — +, q + ~ 2/c and g_ = 0, yielding for both 
i/ = ± the analytical result 
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with the Heaviside function 9. The Eliashberg functions 
for the other two phonon modes have to be computed 
numerically. Together with Eq. ( 13 1 they are shown 
in Fig. [l] Numerically, after a rescaling we find al- 
most universal behavior in the sense that the functions 
(kF/k) 2 F^ u \kSl/kF) are essentially independent of k. 
The O — j- behavior is dominated by the Rayleigh mode 
with F(Si) ~ SI 2 , but at higher energy scales (in partic- 
ular outside the regime shown in Fig. [TJ , the two other 
modes are much more important. 



The resulting quasi-particle decay rate T^ (T) then fol- 
lows from Eq. ( 10 ) and is shown in Fig. [2} The decay rate 
is dominated by the L mode except for very low energy 
scales, i.e., when the particle is near the Fermi surface, 
k ~ hp, and temperature is low, T < Tbg- For high 



temperatures, however, Eq. ( 10 1 predicts a characteristic 
r T law, which allows to identify electron-phonon scat- 
tering processes in practice. For k = kp and T -C Tbg, 
the decay rate is dominated by the R mode, and we ob- 
tain 
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with C(3) ~ 1.202 [26]. This T 3 law and the crossover 
to the linear T dependence for T 3> Tbg are shown in 
Fig. [2j We note that away from the Fermi surface, the 
T = decay rate stays finite and scales as ~ \k — kp\ 3 
for k — » kp. 

Resistivity. — Next we compute the phonon contribu- 
tion to the resistivity, p, using a quasiclassical Boltz- 
mann transport theory as employed recently for graphene 
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with the density of states D(E) = \E\/(2irvp). This ap- 
proach is valid for |/x|(r) 3> 1, which is equivalent to 
Gqp -C 1 with the conductance quantum Gq = e 2 /h. 
The inverse of the energy-dependent electron-phonon 
transport scattering time r(eks) follows from Fermi's 
golden rule as a sum over independent phonon mode 
(A) contributions. The result can again be expressed 
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FIG. 3. (Color online) Temperature dependence of the 
phonon contribution to the surface resistivity. Note the 
double-logarithmic scales. Individual contributions of each 
phonon mode are also shown. 



using a "transport" Eliashberg function (fi) given 

by Eq. (11) with ui = e ks and an additional factor 



(1 — cos9 kq ) in the integral, where 6 k , q = k+q — k 
is the angle between k and k + q. After the angular 



integration, we obtain T as in Eq. ( 12 ) but with an ad- 
ditional factor (q 2 — q 2 _)/[2k(q + — fcjf in the integrand. 
The resulting functions are depicted in the inset of Fig. [T] 



After some algebra, we arrive at [2"7] 
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In Fig. [3j we show the full T dependence of the phonon- 
induced resistivity p. For T <C Tbg, the resistivity is 
dominated by the O — ^ behavior of ^(£1). The latter 
comes from the R mode with JF ~ f2 4 , which implies 
p ~ T 5 as T —> 0. The prefactor can be evaluated exactly, 
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where £(5) ~ 1.037 26J. We thus recover the standard 
BG power law, p ~ T 5 , as in bulk 3D metals [35] which 
is here caused by the coupling to the Rayleigh surface 
phonon mode. For T Tbg, on the other hand, we 
find a p ~ T law predominantly due to the L mode. For 
T ~ Tbg j all three phonon modes are important. 

Conclusions. — We have formulated an analytically 
tractable effective low-energy theory of the surface 
Dirac fermion state in a strong TI with deformation- 
potential coupling to acoustic phonons. The influence of 
phonons could be observed as characteristic temperature- 
dependent decay rate T of quasi-particles in ARPES, or 
from their T-dependent contribution to the surface resis- 
tivity. The phonon-mediatcd effective interaction among 
surface fermions can also be attractive at low frequencies, 
possibly allowing for superconducting correlations; how- 
ever, this topic as well as studies of the electron-induced 
modification of phonon properties in this system or the 
physics near the Dirac point (hp = 0) are left for future 
work. We hope that our predictions will soon be tested 
experimentally. 
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